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ABSTRACT 

The neutral hydrogen (H i) and the ionized helium (He n) absorption in the spec- 
tra of quasars are unique probes of structure in the early universe. We present Far- 
Ultraviolet Spectroscopic Explorer observations of the line of sight to the quasar HE2347- 
4342 in the 1000-1187 Angstrom band at a resolving power of 15,000. We resolve the 
He II Lyo absorption as a discrete forest of absorption lines in the redshift range 2.3 
to 2.7. About 50 percent of these features have H I counterparts with column densities 
Njji > 10 12,3 per square centimeter that account for most of the observed opacity in 
He II Lyo. The He n to H I column density ratio ranges from 1 to >1000 with an 
average of ~80. Ratios of <100 are consistent with photoionization of the absorbing 
gas by a hard ionizing spectrum resulting from the integrated light of quasars, but ra- 
tios of >100 in many locations indicate additional contributions from starburst galaxies 
or heavily filtered quasar radiation. The presence of He n Lya absorbers with no H I 
counterparts indicates that structure is present even in low-density regions, consistent 
with theoretical predictions of structure formation through gravitational instability. 
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The intergalactic medium (IGM) is the 
gaseous reservoir that provides the raw ma- 
terial for the galaxies that dominate our view 
of the visible universe. By observing distant 
bright objects such as quasars, we can explore 
the IGM by examining the absorption features 
it imprints on the transmitted light. These 
absorption features trace structure in the uni- 
verse at epochs intermediate between the ear- 
liest density fluctuations seen in the cosmic 
background radiation and the distribution of 
galaxies visible today. The distribution of ab- 
sorption features according to redshift (z) and 
the column densities of gaseous material in dif- 
ferent ions reveal the structure of the IGM and 
its density and ionization state. From the ion- 
ization state of the gaseous species, we can also 
infer the processes responsible for ionizing the 
gas, e.g., radiation from quasars in the early 
universe, or from early bursts of star forma- 
tion. 

The lack of smooth Lya absorption by H I 
in quasar spectra (_/) led to the conclusion 
that any diffusely distributed gas must be too 
highly ionized to be visible. Discrete absorp- 
tion features were so plentiful, however, that 
they were dubbed the "Lyman a forest" (2). 
Because ionized helium (He + ) is more difficult 
to ionize than hydrogen (ionization potential of 
54.4 eV versus 13.6 eV) and recombines faster, 
its abundance in the IGM is expected to be 
higher. Searches for a diffuse component of 
the IGM therefore concentrated on absorption 
from He n Lya A304 A (3). This absorption 
was first observed using the Hubble Space Tele- 
scope (HST) (4) along the line of sight to the 
quasar Q0302-003 (z = 3.29). The observation 
showed the IGM to be essentially opaque at 
redshifts of z ~ 3, and it seemed possible that 
the material in the H I Lya forest was suffi- 
cient to account for the He n opacity ( 5) . How- 
ever, because no discrete He n features could 
be observed, this required several assumptions. 
Measurement (6) of the He + opacity at lower z, 



using the Hopkins Ultraviolet Telescope (HUT), 
showed a translucent medium with He n opti- 
cal depth r~latz~ 2.4. The HUT observa- 
tions were of insufficient spectral resolution to 
detect discrete He n absorption features. How- 
ever, models of the sightline towards the quasar 
HS1700+64 = 2.72) required a substantial 
contribution to the He n opacity from He + that 
was more smoothly distributed than the known 
H I Lya absorbers (6, 7) and improved HST 
observations of Q0302-003 also required more 
He II opacity than that predicted by the H I 
absorption (8). Resolution of these inconsis- 
tencies requires observations at higher spectral 
resolution. 

Recent theoretical studies view the high- 
redshift IGM as a tracer of cosmic structure 
formation by gravitational instability. In such 
a scenario, diffusely distributed baryonic mate- 
rial (the protons and neutrons that constitute 
ordinary matter) responds to the gravitational 
influence of the underlying dark matter. The 
most overdense regions collapse first to form 
the earliest galaxies. This leaves the remain- 
ing gas as the IGM. Rather than being a uni- 
form medium filling the space between galaxies, 
the IGM itself should show structur on scales 
larger than that of individual galaxies. Both 
numerical calculations (9-16) and analytic the- 
ory (17-21) link the evolving structure of the 
IGM to the H I and He n absorption visible in 
the spectra of high redshift quasars. 

Given the unsettled observational issues and 
the need to test the theoretical expectations, 
we planned Far Ultraviolet Explorer (FUSE) 
observations of the He n absorption toward the 
bright (visual magnitude V = 16.1) z = 2.885 
quasar HE2347-4342 at high spectral resolu- 
tion. Previous observations of HE2347-4342 
at longer ultraviolet (UV) wavelengths showed 
it to be one of the brightest candidates for 
such observations (22, 23). As with the other 
quasars observed with HST, the He II absorp- 
tion in this object is mostly opaque, but there 
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are wavelength intervals of high transmission 
that suggest we can see the beginnings of the 
He + re-ionization in the IGM in the redshift 
range accessible toward this quasar (22). Be- 
cause He II Lya at z = 2.885 is redshifted to 
1181 A and HST is sensitive only down to 1150 
A, which corresponds to z = 2.79, the HST 
data do not cover a large range in redshift. We 
therefore hoped to trace this re-ionization pro- 
cess over a greater range with FUSE because 
its short-wavelength sensitivity would let us 
probe the IGM to much lower redshifts, down 
to z ~ 2.0 in principle (24). 

We observed HE2347-4342 with FUSE in 
two separate campaigns. The first observation, 
from 17 to 27 August 2000, comprised 351,672 
s, of which 192,610 s was during orbital night 

(25) . The second, running from 11 to 21 Octo- 
ber 2000, accumulated 249,717 s, with 183,630 
s during orbital night. For each observing cam- 
paign, HE2347-4342 was centered in the 30 arc 
sec x 30 arc sec apertures. To maintain the 
optical alignment of the four channels during 
these campaigns, we offset FUSE every other 
day to the nearby UV-bright star WD2331-475 
and performed routine adjustments to the mir- 
ror positions. 

At the faint flux levels presented by HE2347- 
4342 (F\ ~ 3 x 10" 15 erg cm' 2 s" 1 A" 1 at 1200 
A), an accurate extraction of a background- 
subtracted, calibrated source spectrum from 
the two-dimensional (2D) recorded data re- 
quired customized processing. Because the 
LiF channels have a throughput that is ~3 
times that of the SiC channels, we discuss 
only the LiF portions of the spectrum. Be- 
cause there are spatial variations in the back- 
ground that affect the extracted flux level in a 
source as faint as HE2347-4342, we turned 
off the standard background subtraction in 
the FUSE data-processing pipeline and ex- 
tracted a "source + background" spectrum 

(26) . We then manually extracted back- 
ground spectra from the geometrically recti- 



fied detector images using regions lying ad- 
jacent to the source spectrum. Dead spots 
and airglow lines from other FUSE apertures 
were masked out of these background spec- 
tra. We computed a linear slope joining the 
two background regions, and scaled the fi- 
nal result for each of the four detector seg- 
ments to match the level at the location of the 
extracted source+background spectrum. We 
subtracted these background spectra from the 
pipeline-extracted source+background spectra, 
and then used the rest of the normal pipeline 
processes to perform the standard wavelength 
and flux calibrations. At all steps we propa- 
gated a 1-cj error array along with the data. 

The multiplicity of detectors in FUSE and 
our two separate observations allowed several 
checks on our reduction process. At each 
wavelength we obtained two separate extracted 
spectra from each observation for a total of 
four independent spectra. All spectral features 
noted in our combined spectrum were visible 
in each independent spectrum. As a check on 
our flux levels, we computed the scatter among 
the fluxes from the 1000 to 1080 A spectra 
(2.7 x 10~ 16 erg cm" 2 s" 1 A" 1 ) and the 1100 to 
1187 A spectra (2.1 x 10" 16 erg cm" 2 s" 1 A" 1 ). 
For both, the scatter is higher than our l-cr 
errors, indicating some residual systematic er- 
rors. These errors are ~5 to 10% of the ex- 
trapolated continuum flux for HE2347-4342, 
and they are indicative of our overall uncer- 
tainty. As a check on our zero levels, we note 
that the interstellar absorption line C n A1036 
is saturated in all FUSE extragalactic spectra. 
All four extracted spectra containing this line 
have zero flux at line center to within the l-cr 
error bars. At wavelengths of >1150 A, pre- 
vious HST observations of HE2347-4342 show 
several regions where there is little or no flux. 
Again, our extracted spectra also show no net 
flux in these regions to within our l-cr errors. 
Given the consistency among these separate 
observations, we combined the separate spec- 
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tra into a single spectrum with a uniform wave- 
length scale and 0.05 A bins (Fig. 1). 

Because the FUSE bandpass stops short 
of wavelengths where the unobscured contin- 
uum of HE2347-4342 is visible (A > 1190 A), 
we obtained low-resolution HST spectra cov- 
ering 1150-3200 A on 21 August and 16 Oc- 
tober 2000 to establish the continuum level. 
These observations each consisted of a 1060 
s exposure using the Space Telescope Imag- 
ing Spectrograph (STIS) with grating G140L 
and a 600 s exposure using grating G230L. 
HE2347-4342 was slightly fainter (by 7%) in 
the October observation, but otherwise the 
spectra were identical. We scaled the Octo- 
ber observation up to the August flux levels, 
and fitted a simple power law, F\ = 3.31 x 
10- 15 (A/1000 A)" 2 - 40 erg cm~ 2 s" 1 A -1 , with 
an extinction correction to spectral regions free 
of galactic absorption lines. We used a mean 
galactic extinction curve with ratio of selective 
to total extinction Ry = 3.1 (27) and a color 
excess E(B — V) = 0.014 (28), where B is the 
blue-band magnitude. 

The spectrum and extrapolated continuum 
(Fig. 1) shows that for z > 2.72 (A > 1130 
A), the IGM is relatively opaque, with occa- 
sional narrow windows that are nearly trans- 
parent (22, 23). Atz = 2.72, the IGM becomes 
more transmissive, and the opacity systemati- 
cally drops at lower z and shorter wavelengths. 
In this translucent region, it is also apparent 
that discrete absorption features account for 
most, if not all, of the opacity. We have re- 
solved the He n absorption into a He n Lya 
forest, analogous to the H I Lya forest. 

The quantitative character of the evolu- 
tion in opacity with redshift is consistent 
with our expectations. The mean opacity 
over the redshift interval 2.3 to 2.7 is th c ii = 
0.91 ±0.01 (29), similar to the sightline towards 
HS 1700+64 (6). We compare the measured 
opacities to a model (21) of discrete clouds 
photoionized by the general quasar popula- 



tion, valid for epochs after He n re-ionization 
is complete. In the model, ionization fluctua- 
tions on scales of ~4000 km s _1 (~15 A, the 
size of our largest bins) lead to the predicted 
range of opacities (Fig. 2). For z < 2.72, where 
He II reionization is potentially complete, the 
absolute level of the observed opacity, its trend 
to lower values at lower redshift, and the fluc- 
tuations about the mean are all consistent with 
this model. 

The extrapolated STIS continuum is consis- 
tent with the peak fluxes in the FUSE spectrum 
(Fig. 1). Positive residuals from points above 
the extrapolated continuum have a Gaussian 
distribution consistent with the la error bars 
on our data points, indicating that there are 
absorption-free windows where the IGM is 
transparent in He II Lya. For models in which 
the physical density of smoothly distributed 
He II evolves as nneii o< (l + z) a , and assuming 
a deceleration parameter A = 0, and density 
parameter Q, = 1, the opacity in smoothly dis- 
tributed He II Lya varies with z as T smoot h(-z) = 
T S mooth(^ = 2.885) [(1 + z)/(l + 2.885)]"- 1 - 5 for 
z < 2.885 (1). By requiring that the distri- 
bution of the ratios of positive residuals to 
the 1-0" errors be consistent with a Gaussian 
of mean and dispersion 1, we set a 1-a up- 
per limit on T smoo th(2 = 2.885) by permit- 
ting T smoot h(z = 2.885), the power-law in- 
dex, and its normalization to vary freely un- 
til x 2 increases by 1. This yields an upper 
limit of T smoo th(^ = 2.885) < 0.11 to 0.12 for 
a = to 6. 

Because the FUSE data resolve the He II 
Lya forest, we can make a detailed comparison 
of the individual absorption features to their 
counterparts in the H I Lya forest. To measure 
the individual He n column densities, we fit the 
FUSE spectrum using the program SPECFIT 
(30). Our model incorporates the extrapolated 
continuum, individual He II absorption lines 
with Voigt profiles (assuming Doppler widths 
of £>Heii = bui), and foreground interstellar ab- 
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sorption in the Milky Way (31). For each He n 
line, we permit the column density to vary, but 
we fix z and the line width (&Heii = &hi) to 
have the same values as those in the H I line 
list from a Keck observation of HE2347-4342 
with a resolution of 8 km s _1 (32). The H I 
lines identified in the Keck spectrum (184 in 
total) all have He n counterparts (Fig. 3). Al- 
though these lines account for most of the ob- 
served opacity in He n Lya, they account for 
only ~50% of the observed spectral features. 
To model the remaining features, we inserted 
additional lines whose redshifts were permit- 
ted to vary freely and whose widths were fixed 
at b = 27 km s" 1 , the mean of the distribu- 
tion from the H I line list. Over the 1000 
to 1130 A wavelength range, 179 additional 
He II Lya absorption features were added to 
our fit. At the 3-a confidence limit, we are 
sensitive to He n lines with a limiting column 
density of log N (Hell) > 10 12 - 8 cm -2 . For com- 
parison, there are 72 H I Lya features with 
Nhi < 10 13 cm -2 in the same corresponding 
range in z in the Keck spectrum. 

The added features are consistent with an 
extrapolation of the number of H I lines per 
unit column density of the form /(Nhi) °t 
N H }- 5 (21, 33, 34) down to column densities 
of ~ 10 11 cm" 2 . This is potentially a substan- 
tial amount of material, but, given the power- 
law form of the column density distribution, it 
does not dominate the total mass in the Lya 
forest (19). Because the physical sizes of the 
absorbers and the ionization corrections are un- 
certain, we refrain from estimating the amount 
of mass that this population of absorbers repre- 
sents. However, in the context of hydrodynam- 
ical models of the IGM (14, 19), we note that 
the low-column-density extension of the Lya 
forest contains ~10% of the baryons at z = 2 
to 3. Our detection of a large number of He II 
absorption features with no H I counterparts is 
consistent with the predictions of these models. 

With measured ratios of He n to H I col- 



umn densities, rj = N (Hell) /N (HI), we can in- 
fer the shape of the ionizing spectrum illumi- 
nating the absorbing gas. Because our data 
include both measured values and lower limits 
(Fig. 4) (for He n Lya features with no H I 
Lya counterparts), we use the Kaplan- Meier 
estimator for censored data (35) to derive a 
mean value (log rj) = 1.89 ± 0.04 for the full 
data set. Models of an IGM photoionized by 
the integrated light from quasars propagated 
through the IGM (21, 36) predict values of 
i] = 30 to 100 for quasar spectra with spectral 
indices a q = 1.5 to 2.1 (f u oc v~ aq ). Intrin- 
sic spectral indices for quasars (as measured 
down to a rest wavelength of ~350 A) lie in 
this range (37). Thus, most of the observed ab- 
sorption features are consistent with photoion- 
ization by quasar radiation. However, about 
40% of the absorbers have r] > 100, indicating 
that localized regions of the IGM are photoion- 
ized by softer spectra. This might be heav- 
ily filtered quasar radiation (in regions where 
the re- ionization of He n is not yet complete), 
or stellar radiation from ongoing star forma- 
tion (21) 

There are fluctuations in r\ on the scale of in- 
dividual spectral features (Fig. 4) which implies 
that the ionizing radiation field is not uniform. 
Local radiation sources (which may be either 
quasars or star-forming galaxies) probably af- 
fect T) and the local opacity of the IGM. This 
characteristic is most evident at high z in the 
low-opacity "windows" first seen in EST spec- 
tra of HE2347-4342 (22, 23), but it persists to 
lower redshift where the IGM is less opaque. 

By resolving the He II Lya absorption in 
the IGM using FUSE, we are able to see struc- 
ture in the universe that extends into the low- 
est density regions. The presence of He n 
absorption with no H I counterparts is con- 
sistent with hydrodynamical models that pre- 
dict density fluctuations due to gravitational 
instabilities on all scales, from the high den- 
sity peaks that form galaxies to the distribu- 
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tion of gas in low-density voids {9-16, 19, 20). 
Observations of He n Lya absorption is the 
best method to study these regions, and fu- 
ture observations of additional bright quasars 
with FUSE should provide essential informa- 
tion on the cosmic variance in the structures 
we see along different lines of sight. 
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Fig. 1. — FUSE spectrum of HE2347-4342 (binned to 0.05 A pixels) is shown in black. Bins at the 
peak fluxes have a signal-to-noise ratio of ~7. A portion of the contemporaneous STIS spectrum is 
shown in green. The red line is the extrapolation of the power law plus extinction of E(B-V)=0.014 
fitted to the STIS spectrum. The position of He II Lya at the redshift of HE2347-4342 is marked. 
Gaps in the FUSE spectrum due to excised terrestrial airglow lines are marked with green 0. The 
broad gap from 1072 to 1089 A is due to gaps between the FUSE detector segments. 
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Fig. 2. — The He II opacity THeii m coarse 5 to 15 A bins is shown as a function of redshift. The 
solid curve is the theoretical prediction for the opacity assuming line blanketing only (21). The 
dashed curves indicate the 1 a extent of variations in the opacity anticipated due to ionization 
fluctuations. At redshifts below z = 2.72, where the reionization of He II in the IGM appears to be 
complete, the model matches well the observed opacity, its trend to lower values at lower redshift, 
and the fluctuations about the mean. 
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Fig. 3. — (A) The black line is a portion of the FUSE spectrum of HE2347-4342. The smooth red 
curve across the top is the extrapolated continuum. The magenta line overlaying the spectral data 
in black is the model described in the text. The area shaded in yellow shows the fraction of the 
He II opacity due to absorption features that correspond to H I absorption lines identified in the 
Keck spectrum (32). The area shaded in green shows the fraction of the opacity due to additional 
He II absorption features that have no H I counterparts in the Keck spectrum. (B) The normalized 
Keck spectrum (32). Note the direct correspondence between the wavelengths of H I lines in the 
Keck spectrum and He II absorption features in the FUSE spectrum. 
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Fig. 4. — The logarithm of the He II to H I column density ratio r/=N(HeII)/N(HI) versus redshift 
for measured absorption features in the FUSE spectrum. The black points with error bars have 
measured column densities for both Hell (from FUSE) and H I (from Keck). The red points are 
lower limits computed using the He II column density measured from the FUSE spectrum and 
an upper limit of 10 12 ' 3 cm -2 for H I absorption lines in the Keck spectrum. A Spearman rank 
correlation coefficient of 0.05 indicates that there is no significant trend with redshift. The mean 
value, (log 77) = 1.89, is consistent with photoionization of the IGM by integrated quasar light {21, 
36). However, one can see significant fluctuations in 77 that indicate that the ionizing radiation field 
is not uniform. Values of rj > 100 indicate regions photoionized by either heavily filtered quasar 
radiation or bursts of star formation. 
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